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Abstract. The optical conductivity of charge carriers coupled to quantum phonons is studied 
in the framework of the one-dimensional spinless Holstein model. For one electron, variational 
diagonalisation yields exact results in the thermodynamic limit, whereas at finite carrier density 
analytical approximations based on previous work on single-particle spectral functions are 
obtained. Particular emphasis is put on deviations from weak-coupling, small-polaron or one- 
electron theories occurring at intermediate coupling and/or finite carrier density. The analytical 
results are in surprisingly good agreement with exact data, and exhibit the characteristic 
polaronic excitations observed in experiments on manganites. 



PACS numbers: 78.20.Bh, 71.10.-w, 71.27.+a, 75.47.Lx, 78.20.-e 



1. Introduction 

In recent years, there has been a lot of controversy about the possible interpretation of optical 
data on strongly correlated electron systems such as manganites, cuprates and nickelates in 
terms of polaronic charge carriers [1,2]. It turns out that simple one-electron or perturbative 
weak/strong-coupling theories can usually not explain more than one type of experiment for 
a given set of (realistic) parameters, and a general discussion of this issue for the case of the 
manganites has been given in [3-5]. Owing to these discrepancies, and also in the light of 
recent progress in the understanding of many-polaron systems [6-8], it is highly desirable to 
revisit this problem using both analytical and numerical many-body methods. 

Ignoring the necessity of taking into account cooperative Jahn-Teller effects to describe 
orbital effects in manganites [4], an important issue in the framework of one-band models 
concerns the validity of two complementary types of electron-phonon interaction, namely the 
Holstein [9] and the Frohlich model [10] for screened, local coupling and unscreened, long- 
range electron-phonon interaction, respectively. Motivated by recent experimental results 
and the absence of a unified theory for, e.g., ferromagnetic manganites [5], we analyse the 
signatures of polaronic excitations in the optical conductivity, both in the low-density limit 
and at finite carrier density, in the framework of the one-dimensional spinless Holstein model. 

Polaronic materials such as the manganites generally require a theory valid for finite 
carrier density and all coupling regimes [4,5]. Whereas a large number of works have 
reported on approximate findings for the one-electron limit (the polaron problem) of the 
Holstein model, here we present variational diagonalisation results for the thermodynamic 
limit, thereby providing an exact (numerical) solution. Exact data for finite clusters, also at 
finite temperature, can be found in [11]. 
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The many-electron case has been studied by approximate methods in the past [3, 12- 
14]. Moreover, exact results on finite clusters are available [15,16]. In order to further 
improve the present understanding, we examine the optical conductivity in the framework of 
electronic spectral functions deduced by perturbative and variational methods valid at finite 
carrier density [17, 18]. The important advantage of such calculations is the possibility to 
relate the observed contributions to the optical response to specific processes (transitions) in 
the analytical formulas, and to understand the relevance of these contributions in dependence 
on the model parameters. Special attention will be paid to the role played by the coherent 
(respectively incoherent) parts of the spectral functions. We compare the analytical results 
to exact numerical data, and discuss the differences between the weak-coupling (WC), 
intermediate-coupling (IC) and strong-coupling (SC) regimes, as well as the influence of 
carrier density and phonon frequency. Besides, we attempt to make a qualitative connection 
of our findings to recent experiments on manganites [5], as well as to theoretical work on 
the many-electron case based on the Holstein model [3, 16] and the Frohlich model [19, 20]. 
Theoretical work on manganites is discussed in [3-5], and a review of polaron theories can be 
found in [21]. 

This paper is organised as follows. In section 2, we introduce the spinless Holstein 
model. In section 3, we outline the derivation of the general expression for the optical 
conductivity involving the electronic spectral functions deduced in [17]. Numerical and 
analytical results are discussed in section 4, and section 5 contains our conclusions. 



We consider the one-dimensional (ID) spinless Holstein model, describing fermions coupled 
to dispersionless optical phonons. It provides a general framework to study polaron physics 
[22], many-polaron effects [8] and quantum phase transitions [23-26], but is simple enough to 
permit reliable investigations by analytical and numerical methods. A connection to the more 
general Holstein-Hubbard model (see, e.g., [15]) can be made in the limit of large Hubbard- 
U in the latter. Moreover, the Holstein double-exchange model for the manganites [3,4,27] 
reduces to the present model for infinite Hund's rule coupling and a ferromagnetic state. 
Following [17], we write the Hamiltonian in the general form 



where c] (c i ) and b\ (b i ) create (annihilate) a spinless fermion respectively a phonon of energy 
u>o (H = 1) at site i. The strength of the electron-phonon interaction is specified by the 
dimensionless coupling constants A = Ep /2t (adiabatic regime, oj /i « 1) and g 2 — Ep /uq 
(anti-adiabatic regime, wo/i » 1), with Ep denoting the atomic-limit [CV, = for i ^ j in 
equation (1)] polaron binding energy. 

The definitions of r\ and the coefficients dj will depend on the type of approximation 
used and hence on the parameter regime. In the WC case, in which we start with the original, 
untransformed Holstein Hamiltonian, we have 



where fj, denotes the chemical potential. In contrast, in the SC regime, we use the Hamiltonian 
after the Lang-Firsov transformation [28] with 



2. Model 




(1) 



»? = -M, Cu =0wo(&j + &i), C(ij)=t, 



(2) 



v = -Ep-f,, c u = o, c (ij) = te-»< 6 !- 6 «- 6 H^ . 



(3) 
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The non-interacting (polaron) half-bandwidth in the SC regime is W = Wc~ 9 , with the ID 
free-fermion half-bandwidth W = 2t. Finally, the treatment of the IC case is based on a 
modified Lang-Firsov transformation [17] defined by the unitary operator 

JJ = e E» ff(7ctc i +7)(i>t-6 i ) ^ ^ 

which leads to a Hamiltonian of the form (1) with 

r] = — fi — E P [7(2 — 7) + 2 7 (1 - 7)] , 

C u = . 9 co(l - l){b\ + K) , C (ij) = tc--^- b > - b l +b ^ . (5) 

The parameters 7 and 7 depend on the carrier concentration n, and on the variational 
parameter R characterising the delocalisation of the charge carrier from the center of the 
associated lattice distortion. Explicitly, we have 

7 = 2nc~ 1/R tanh ^ , 7 = tanh ^-7, (6) 

and R is then defined by the position of the minimum of the total energy per site in the first 
(Hartree) approximation [17]. 



3. Theory 

In this section we present our analytical results for the optical conductivity. We first 
derive a general expression which allows us to calculate the real part of the optical 
conductivity from the momentum and energy dependent one-electron spectral function 
A(k 1 10). Approximations for the spectral function in different coupling regimes have been 
obtained in [17], and the basic formulas necessary for the computation of the conductivity are 
compiled in Appendix A. The calculations presented in [17] are based on the self-consistent 
treatment of the self-energy equations in second-order perturbation theory. To describe the 
crossover from the SC to the WC regime, we use the variational procedure outlined above. 

Our approximation neglects vertex corrections due to direct fermion-fermion interaction, 
which is also absent from the spinless Hamiltonian (1). Nevertheless, fermion-phonon and — 
at finite carrier density — phonon-mediated, retarded fermion-fermion interaction effects enter 
via the spectral functions [17]. For simplicity, we exclude from our discussion the quantum 
phase transition from a Luttinger liquid to an insulating Peierls phase at half filling [23, 25]. 



3.1. General expression for the optical conductivity 

The current density operator for the ID spinless Holstein model reads 

^ =i T7 H id' -9)cl>c g , (7) 
(g',g) 

where g', g number the lattice sites, and a = \g' — g\ respectively V = Na are the lattice 
constant and the volume. 

According to linear response theory [29], the frequency dependent complex conductivity 
ct{lu) is determined by the Green function of the current density operators. In particular, 

Recr(w) = --ImG?(w) , (8) 

u> J 

where 

/oo 
d(t - t')Gf(t - ty^-O (9) 
-00 
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with the retarded Green function 

G f(t-t') = -i([j{t),m})o(t-t'). do) 

In the notation of [30], the related Matsubara Green function takes the form 



H (rm - m 2 )(m' 1 - m' 2 

(mi,m2) (m'^mj) 

x^rd^n^^njct {r[)c M)) 



mi 

2 



(J) H H ("»i-m2)(mi-n4; 

(mi, 7712) (m^,m 2 ) 

x (T r c m2 (r 2 )c m ,(^)ct (rOc^Cn)) 



,(H) 

t 2 =t 1 -,t! ! =t[- 

where rrii, m\ are integers labeling the sites of the ID lattice. The two-particle Green function 

G(2,2';l,l') = (T TCro2 (75)^, (7^)4, (7^)4, (n)> (12) 

for mutually independent fermions can be written in terms of one-particle Green functions 
[30]. Assuming (j) = in the absence of an electric field, the part of the two-particle Green 
function relevant for the evaluation of equation (11) reads [30] 

G(2,2';l,l')= -G(2,1')G(2',1) 

= - (T T c m2 (r 2 )cl,(r[))(T T c m , (r'^in)) . (13) 

We use the Fourier transformation of the one-particle Green functions 

G(2,l') = ^^e ife(m2 -" i 'i ) i^e- i ^ (T2 -^ ) G(fc,i^) (14) 

k u>„ 

with lo v = [2v + 1)tt//3, k = 2nm/N (likewise for G(2', 1)), and the spectral representation 

G{k,iu v ) = r dcj' A(KLU '] , (15) 

to obtain 

rP 

Gj{iw n ) = / d(n - riy^-^Gjin - t[) , (16) 
Jo 

where uj n — 2nir/fi. Using equations (12)— (15) to express the rhs of equation (11), the 
Fourier transformation (16) is evaluated by carrying out the summations over m,i, m' i , and the 
integration and subsequent summation over the Matsubara frequencies lo u [30]. We find 

= f^y E( sinfc ) 2 r ^ r ^^^^ ^y^it? ■ a?) 

Here f(x) = [exp(j3x) + is the Fermi function. Using Dirac's identity 

- ; 1 ; - - n+ = — -, --iTr^ + w'-e'), (18) 

iv + iv' — e + i0+ iv + iv' — e 

the analytical continuation \u n 1— » iv + iS gives in the limit 5 — > + 

ImGf(w) = -it ^(sinfc) 2 / dw' / de'A(fc, e') 

x[/( W ')-/(e')]<^W-e'). (19) 
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In particular, for uj > and T — > 0, and introducing the Heaviside step function 9(x), 

(Opf \ 2 
—J ^{sinkf J dw'ACjfe.w'J^fe.w' + wJ^w' + w). (20) 

The real part of a(u>) for w > 0, given by 

a^(oj) = ^°-L V(sinfc) 2 / dw^fc.w'^Cfe.w + wO^w+w'), (21) 

is determined by the overlap of the electronic spectral functions for energies uj' < (relative 
to /li) and uj' + u) > 0, respectively. Here we have defined ct = a7r(et) 2 . 

The neglect of vertex corrections in the current-current response function based on 
equation (13) leads to the factorisation into spectral functions in equation (21), which is 
equivalent to the exact D = oo result commonly used in dynamical mean-field theory (DMFT) 
[31]. However, in contrast to the local approximation of DMFT, the spectral functions used 
here have a non-trivial momentum dependence [17]. 



3.2. Limiting cases 

Weak coupling The explicit form of the WC electronic spectral function, calculated 
according to equations (A. 1 )-( A.4), was determined in [ 1 7] . Taking into account the condition 
\uj\ < uj n (\ui\ > uj ) for the coherent (incoherent) part, and the explicit result (A. 2) for the 
coherent part A c (k, uj), we find for the optical conductivity 



<7 res (u/) = — / dfc(sinfc) 





z k 0(-(E k - n))6(w +E k - n)6{E k ~n + u- tu )A ic (k, E k -/j, + uj) 
+z k 6(E k - h)6(uj - (E k - n))9{ui - (E k - p) - oj )A [c (k, E k -fi-uj) 



r° 

+ / duj'A lc (k, iu')A lc (k, u + J)6{-u' - lu )0(uj +uj' - oj q ) 

J — OO 



(22) 



For the discussion of equation (22) in section 4, we have to recall [17] that the WC 
approximation of the incoherent one-electron spectral function A lc (k,u') is non-zero only 
for uj' G (—ujo — W — /i, —ujq) or u>' £ (luq, ujq + W — /i). The last term in equation (22), 
corresponding to transitions between these frequency intervals of the spectrum, is denoted in 
figure 1 as D. 



Intermediate coupling Within the variational approach, the electronic spectral function is 
given by equation (A.5) with g replaced by ~fg and W — We' 1 9 . The explicit expression 
for T,(k,uj) determining the polaronic spectral function can be found in [17]. Using the 
definitions (A.6)-(A.9), we can write 

r0 

&uj'6{uj + uj')A c (k, uj')A c (k, uj' + uj) = A + B + C + D . (23) 



/ 

j — ( 



In general all terms of equation (23) contribute to er rcg (u;). We find 

4cr n 



TT 2 U) ^ S\ 

s>l 



dfc(sinfc) z k 



o 



\d k , E k 



r0(E k + r,)8(uj - (E k + v)M-E kl - V )9(lo + E kl + r,) 
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Zk- 2 



\d k2 E, 



k 2 \ 



-E k - r,)6{uj + E k + V )d(E k2 + vMloo - {E k2 + r,)) 



2a 



+— e 



E 



2 s (73) 2s 



si 



-(1 



r 7T 

-2 1 -) / 
Jo 



dkz k - 



z k > 



d k 'E k / 



x6(-E k - V )6(u> + E k + v )9{E k , + v,)6{w a - {E k , + r,)) , 



(24) 



where 

z k 



\d k E k \ 



W sin k J 1 - - cos k 



S>1 



w 



w 



d£ 



i - (a/wy 



E k - sluq - £ E k + suj - £ 



xV 



/w 
-w n 



d£ 



0(Z + V ) 



H — 7(1 — "f)E-p sin A; 

7T 



E k — loq — £ E k + uj — £, 



(25) 



Here fci,/c2 G (0,7r) are determined by the conditions E kl = E k — (uj — sujo) and 
E k2 = E k + (uj — sojq), whereas k! > kp is given as the solution of Ey = E k + (uj — suo). 
For the other contributions to equation (23), we have 

B = c-^ 2 z k Af(k,iu + E k + V )6(oj + E k + ri)6{-E k - V )6(oj + E k + rj) 

(jg) 2s 1 



+c 



-m 2 



E 



S\ 7T 



dk'z k 'A^(k, uj + E k ' + V - slo ) 



C 



x 6(u) + E k '+t]- su) )9{-E k , - r))0(E k > + t] + uj ), 
-^ 2 z k A^(k, E k + V - uj)6(w -E k - r{)6{E k + 7?)0(w o -E k -rj) 

+C -(7 S ) 2 y-Mli f dk'z k ,Ai c (k,E k , +TJ + SCO0-U) 

x9(uj - E k , - 7) - suj )9(E k , + r))6{u) - E k , - 77) , 



(26) 



and 



D = 



f 

J —1 



duj'A^(k, uj')A^(k, uj + u>')6{u) + uj') . 



(27) 



(28) 



so that 



[g^(uj)] b+c+ d = — r dfc(sin k) 2 (B + C + D) . (29) 

Jo 

The parts of the electronic spectral function denoted by Af(k, uj'), defined in equation (A.8), 
are determined by the incoherent parts of the polaronic spectral function. 



Strong coupling For g 2 3> 1, small polarons are the correct fermionic quasiparticles, and 
the calculations corresponding to equations (A.1)-(A.4) are done using the small -polaron 
Hamiltonian defined by the coefficients (3) (see [17] for explicit results). In this regime, the 
zeroth-order small-polaron approximation — with the spectrum consisting only of coherent 
states in the polaron band — is applicable and the polaronic spectral weights z k « 1 [17]. 
Consequently, it is sufficient to keep only the term A, and SC expressions corresponding to 
equations (24) and (25) can be obtained by setting 7=1. 
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Figure 1. Illustration of the transitions contributing to equation (22). 
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Figure 2. Transitions contributing to equations (24), (26) and (27). Equation (28) corresponds 
to figure 1(D) with — fi replaced by i). 



The rather complicated explicit formulas for the optical conductivity presented above are 
visualised in figures 1 and 2 for the WC and the SC/IC cases, respectively. In particular, these 
pictures show that non-zero values of <t ics (lj) are restricted to uj > ui n . The existence of such 
a threshold for the optical conductivity at T = was already recognised in [20] on the basis 
of WC calculations in the framework of the continuous polaron model. 



4. Results 



We now come to a discussion of numerical results obtained from the analytical expressions 
derived in section 3, which will be compared to exact data for one electron. It is instructive to 
first consider the kinetic energy and the Drude weight, related to cr res (w) by the f-sum rule. 
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Figure 3. Numerical results for the kinetic energy E^ n ( ), the Drude weight T> ( ), 

and the integrated spectral weight S rcg ( , equation (33)), for one electron as a function 

of electron-phonon coupling strength. Here (a) wn/t = 0.4 and (b) uio/t = 4. 



4.1. Integral quantities 

In a tight-binding system with the kinetic energy operator 

and the current density operator (7), it follows that [32, 33] 

— e 2 a£kin = — / da)cr(cj) = — / dwRecr(w), (31) 
n J-oo n Jo 

by symmetry of u(lo). Here, £ki n = (T) /N represents the kinetic energy per lattice site. 

For uj > and T = 0, the electron-phonon interaction implied in the Holstein model 
considered here gives a non-zero Re a{to) only if lo > un- At to = 0, Re <j(u) has a singularity 
to be deduced in Appendix B. Writing 

Re a(uj) = VS(uj) + a lcg (w) , (32) 

where V is the so-called Drude weight, and defining 5 rcg = 5 rcg (oo) with 

S ics (lu) = f duj'a ICS (uj') , (33) 
Jo 

the f-sum rule (31) takes the form 

-e 2 ae ki „ = -P+-5 rcg . (34) 

The kinetic energy may be determined independent of the rhs of equation (34) from the 
electronic spectral function as [34] 

1 f a 

£ kin = ^X!/ &-ue k A(k,u), e k = -2tcosk. (35) 

k J -°° 

In the case of one electron in a lattice with N sites (i.e., with carrier concentration 
n = 1/N), £kin is equal to the kinetic energy E^in of one electron divided by N. To 
obtain non-zero results, both sides of equation (34) are multiplied by N before taking the 
limit N — > oo. Accordingly, using the spectral representation at T = 0, we have for u> > 

aieg{uj)= a^ ^ K^IjI^)! 2 ^-^-^)] (36) 
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Table 1. Analytical results for the kinetic energy — Skin/* f° r wo/t = 0.4, QMC data for 
N = 16, fit = 10, At = 0.05 [7], and ED data for N = 10 (*-not fully converged with 
respect to the number of phonons M = 25) [16]. 



Parameters 


Analytical 


QMC 


ED 


Ep/t = 0.1, n = 0.1 
Ep/t = 0.1, n = 0.3 
£ P /i = 0.1, n = 0.4 


0.202 (WC) 0.211 (IC) 
0.509 (WC) 0.518 (IC) 
0.600 (WC), 0.609 (IC) 


0.190 
0.504 
0.596 


0.1986 
0.5208 
0.5824 


.Ep/t = 2.0, n = 0.1 
Ep/i = 2.0, n = 0.3 
£p/i = 2.0, n = 0.4 


0.155 (IC) 
0.382 (IC) 
0.397 (IC) 


0.136±0.002 
0.36 ±0.01 
0.43 ±0.02 


0.1517 
0.4070 
0.4531 


Ep/t = 4.0, n = 0.4 


0.13 (SC) 


0.19±0.01 


0.2318* 



with the current operator J = -iet J2i( c l c i+i ~ c l+i c i) f° r a single electron. Here \ip m ) 
denotes the m-th excited state with energy E m . In comparison to equations (7), (8) a factor 
l/N is omitted in equation (36). 

According to Kohn's formula, the Drude weight of the Holstein model with one electron 
can be calculated from the effective mass m* as (Ek is the polaron band dispersion) [35] 

V = e\a 8 2 k E k \ =e\a^. (37) 

The evaluation of V via the f-sum rule, and the determination of m* are independent 
computations, so that equation (37) may be used to validate the numerics. 

4.1.1. One electron We adopt the basis construction of [36], which allows for the calculation 
of results that are variational in the thermodynamic limit, and find that the f-sum rule is 
fulfilled to at least six digits for the parameters shown. For spectral properties a Chebyshev 
expansion method is used [37]. 

The dependence of i^kin, S ICS and V on the electron-phonon coupling strength reflects 
the well-known crossover from a large polaron at weak coupling to a small polaron at 
strong coupling [38]. As expected, the results in figure 3 reveal a significant dependence 
on the adiabaticity ratio Uo/t. The adiabatic regime (figure 3(a)) is characterised by rather 
pronounced decrease of i^kin and V (increase of m*) in the vicinity of the point A = 1 near 
which the crossover occurs, which is compensated by an increase of S lcg due to enhanced 
incoherent scattering. In contrast, in the non-adiabatic regime (figure 3(b)), these changes 
occur over a much larger range of the relevant coupling constant g 2 . 

4.1.2. Many electrons Equation (35) turns out to be numerically problematic in the SC 
regime, and £kin is instead calculated from equation (B.17). An approximation for V is given 
by equation (B.10). 

To test the reliability of the spectral functions obtained from the analytical approach we 
compare in table 1 the analytical kinetic energy £ki n for different parameter sets to numerical 
data from quantum Monte Carlo (QMC) simulatons [7] and exact diagonalisation (ED) [16]. 
We restrict ourselves to the adiabatic regime since the analytical approach generally works 
better for Uo/t ^> 1 [17]. The agreement is satisfactory in all cases, with the deviations 
of the analytical results for Ep ft = 2 originating to some degree in the missing spectral 
weight in A(k,oj) (see below). The SC approach slightly underestimates the kinetic energy 
for Ep/t = 4, a coupling which does not fall into the true SC regime (see also figure 5(c)). 
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Figure 4 shows the dependence of £ki n , f , and S res on carrier density for Ep /t = 0.1. 
Whereas the kinetic energy and the Drude weight are almost identical in the adiabatic and 
non-adiabatic regime, the regular part 5 rcg is substantially larger for uj /t < 1 owing to 
enhanced incoherent scattering of carriers by phonons, as is also observed in the one-electron 
case (figure 3(a)). The weak coupling leads to very small S TCg (note the different scale used 
for S Icg ) and hence requires — via the f-sum rule — that — £kin/i ~ tV/ao- A comparison of 
S re& to exact IC results [16] reveals similar deviations as for £ki n - 

The approximations for the Drude weight V derived in Appendix B yield reliable results 
in the WC and SC cases, whereas the dependence of V on Ep for intermediate coupling is 
not properly described. Alternatively, an estimate for T> may be obtained from the sum rule, 
but will be affected by the deviations of £ki n and 5 rcg from exact results. 

4.2. Optical response 

4.2.1. One electron Although the optical conductivity for one electron has been discussed 
before [38], we shall begin with this case as it provides a framework for the finite-density 
case, and another test for the analytical theory. We focus on the adiabatic regime uj /t < 1 
relevant for many polaronic materials. 

Figure 5 shows numerical and analytical results for a rcR {uj) for vo/t = 0.4 and three 
different values of Ep. In order to compare with the one-electron case, we have chosen 
n = 0.1 in the analytical approach, and scaled the results to yield the exact one-electron value 
for 5 reg . The dependence of <j reg (w) and 5 reg on n will be discussed below. 

As anticipated from the analytical results (cf equations (22) and (24)), there exists a 
minimal absorption energy (absorption threshold) u> = luq. Another striking feature is the 
change from an asymmetric spectrum in the large polaron regime to a (more) symmetric 
spectrum in the small-polaron regime also observed experimentally [5]. 

Starting with weak coupling Ep = 0.1 in figure 5(a), we find a pronounced absorption 
signal near uo = uj , followed by a continuous decrease. The analytical results from the IC 
approximation reproduce the main features (notice the good agreement for 5' reg (w)), with the 
sharp low-energy peaks smeared out. As illustrated by the WC results in the inset (in good 
agreement with the IC results), this smearing depends on the number of /c-values used in the 
integration (N = 31 and 801, respectively). The inset also contains analytical results for 
/j, — > Eq, corresponding to the one-electron case [18], for which a non-zero WC result is 
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Figure 5. (colour online) Numerical one-electron results ( , red) for the regular part of 

the optical conductivity cr Tcg (u>) and the integrated spectral weight S rce (uj), and analytical 

results ( , black, scaled to give the same 5 rcg ) for n = 0.1 from the IC approximation 

with (a) R = 2.4, (b) R = 2.1, (c) R = 0.1. Here u) /t = 0.4 and (a) E P /t = 0.1, (b) 
Ep /t = 2, and (c) E-p jt = 4, The inset in (a) compares the WC analytical result for n = 0. 1 

( , black) to that obtained for n — * ( , blue, see text), both rescaled to 5 rog for one 

electron. Panel (c) includes the small-polaron result of [39] ( — ■ — , blue). 



obtained by using the zero-density spectral function in equation (22). Obviously, small peaks 
in <7 reg (a;) are further washed out in this limit. The sharp peaks in the exact numerical results 
are due to the use of a finite number of poles in equation (36). 

Turning to the IC results (Ep /t = 2, close to the critical coupling for the small-polaron 
crossover) depicted in figure 5(b), we find low-energy peaks split off from the high-energy 
part and separated by luq. Besides, the long high-energy tail of both the exact and analytical 
curves is reminiscent of experimental data on manganites [5] and also Ti0 2 [40], and its 
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Figure 6. (colour online) Numerical results for a Teg (uj) for one electron ( , red), 

and analytical results ( , black) for n = 0.1 (scaled to the same S rcg ) from the WC 

approximation. Here u>o/t = 4 and Ep/t = 0.1. Also shown is S Tes (uj). 

occurrence in the IC regime points toward the inadequacy of standard small-polaron theory 
for such materials. 

Finally, for E p /t = 4 (figure 5(c)), we find a maximum close to that of the small- 
polaron result [39] which peaks at 2E p and is included for reference, with the parameter W 
chosen so as to get the same S" cg . For technical reasons, the optical conductivity in the SC 
regime (given analytically by a sum of <5-peaks) is calculated using an artificial broadening 
(see many-electron case). 

The exact numerical results, however, reveal that the parameters chosen do not fall into 
the true small-polaron regime, as indicated by a slight asymmetry, a maximum clearly below 
2E p , and three low-energy peaks (not related to finite-size effects) which diminish with 
increasing E p . The latter originate from transitions between the (dispersive) coherent band 
and dispersionless phonon satellites (cf figure 9(b) of [17]), and are absent in the analytical 
results probably due to an overestimation of band-narrowing in the adiabatic regime. The 
possibility of such overestimation by the Lang-Firsov small-polaron theory in the adiabatic 
case was qualitatively illustrated for the two-site model in [41], and the analytical results of 
figure 5(c) are expected in the SC (small-polaron) regime. 

Figure 6 shows a TCS (uj) at weak coupling in the non-adiabatic regime, taking w Q /t = 4. 
As multi-phonon excitations have small weight due to g 2 <C 1, the absorption is virtually 
restricted to the interval [u> , ui Q + 2W]. Here the analytical WC approach — restricted to one- 
phonon excitations, see equation (22) — yields even better agreement than in figure 5(a). 

Overall, the analytical approach is capable of reproducing the main features of the 
optical conductivity in the low-density regime, including the increase of the frequency range 
for absorption (figure 5). Furthermore, in accordance with figure 3, 5 lreg is small at weak 
coupling, takes on a maximum in the IC regime, and decreases again due to the suppression 
of e kin (via the sum rule) at strong coupling. 

4.2.2. Many electrons There are some issues concerning the analytical results shown 
here which deserve attention. First, we pointed out in [17] that the total spectral weight 
w = J dk J du>A(k, ui) < 1 contained in the spectral function (i.e., its norm) provides a first 
indication for the quality of the analytical approximation for a given set of parameters. 

The IC approach reproduces the WC (SC) approximation in the WC (SC) regime [17]. 
In the vicinity of these limits, it yields results very similar to the WC/SC approximation, with 
w < 1. In the true IC regime (A = E p /2t w 1 or g 2 = E p /uj Q ss 1), the norm w becomes 
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Figure 7. Analytical results for a TC£ (uj) from the WC approximation for n = 0.1 ( ), 

n = 0.3 ( ) and n = 0.4 ( ). Here E P /t = 0.1, (a) wo/t = 0.4 and (b) 

wo /t = 4. The corresponding results for S rcg are given in figure 4. 



as small as 0.7 in the worst case considered (for the parameters of figure 1 1(a)), with w also 
depending on the radius R, whose "optimal" value is determined by minimising the energy 
within the Hartree approach. Despite these deviations, observables such as the kinetic energy 
(see table 1) or the total energy (see figure 13(a) in [17]) are in satisfactory agreement with 
exact results. Moreover, the optical conductivity — calculated from the spectral function — 
is mainly determined by states near the Fermi level, so that missing high-energy incoherent 
processes are irrelevant, and w < 1 does not necessarily lead to poor results. 

For some parameters, the spectral functions obtained from the IC approach exhibits small 
regions in the (k, u) plane where A(k, u) < [17]. To calculate the optical conductivity, we 
have replaced such values by zero. 

Finally, in the SC regime, an artificial "softening" proportional to the energy resolution 
is imposed on the 9 functions in equation (24) which affects the integrated weight 5 reg . 
Nevertheless, a consistent value for S r ° s , and hence the normalisation of cr rcg (w), can be 
obtained using equation (B.16). 

Weak coupling We begin with the density dependence of cr rog (w). Figure 7 shows results 
for Ep/t = 0.1, and we first consider the adiabatic case u> /t = 0.4 in panel (a). The spectra 
exhibit several non-trivial features. Apart from the absorption threshold at oj , we find an 
asymmetric peak exhibiting small but noticeable wiggles which extend over an interval cu 
centred at the maximum, and diminish with increasing n. 

The origin of this feature are transitions of type II, as can be seen by separating the 
contributions I and II to equation (22), and looking at the corresponding spectral function 
A(k, u>) shown in figure 8. 
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Figure 8. (colour online) Coherent (A c , ) and incoherent (A IC , ) parts of the WC 

electronic spectral function for u)o/t = 0.4 and Ep/t = 0.1 [17]. 

Non-zero contributions to the second term in equation (22) exist for oj 6 (ui . , <jJ max ). 
For a given oj from this interval, the range of fc-values contributing to er rcg (w) is defined by the 
condition A lc (k, E k — fi — ui) > 0, with E k — /i > 0. The absorption starts at 0J min = oj q , with 
a single fc-value fc min = fc F . The upper cut-off for the adiabatic case (figure 7(a)) is given as 
w max = + 1 — W— n\, coiTesponding to a single value fc max defined by E k — = uj . In 
the antiadiabatic case (figure 7(b)), for which the entire coherent band lies inside the interval 
(— ui Q , u> ), w max = uj + 2W corresponding to fc max = tt with the highest energy W — \x. 

The shape of the absorption curve in the interval (uJ min , ^ max ) depends on the fc-interval 
contributing to the transition II for fixed oj, and the corresponding values of A lc ( fc , E k — [i — uj ) . 
For the anti-adiabatic case, figure 7(b) shows sharp maxima at the n-dependent positions 
uj 1 = uj q + 1 — W — n\. This feature may be related to the increase of the range of contributing 
fc-values, beginning with the single value k = fc F for u> = u> , up to the interval (fc F , k x ) for 
oj 1 , with k 1 determined by E ki — p, = \ — W — fJ,\. 

For n = 0.1, there is very little weight contained in incoherent states below /i 
(figure 8(a)). In fact, there exists only a single, small and narrow peak with negligible 
dispersion, and transitions from the latter to a coherent state can occur for cj < uj < 2ui Q , i.e., 
for final states between E k — fi = and the upper edge of the coherent band E k — fj, = uj . 
The wiggles, more pronounced for smaller numbers of fc-values, result from transitions to the 
discrete <5-peaks of the coherent spectrum. 

Above ui = 2uj , for n = 0.1, the transitions of type II fall off quickly with increasing 
oj and go to zero for uj > 2uj q + | — W — fi\ w 0.87. In contrast, for larger n, the energy 
interval of incoherent states below the Fermi level | — W — (i\ (and the weight contained 
in the latter) is significantly larger, yielding the cut-offs of u ~ 1.6 and uj rs 2.15. The 
high-frequency optical response is much smoother since the incoherent parts of the spectral 
function are broadened proportional to E p , and 5 I0g (oj) increases in accordance with figure 4. 

Finally, the long high-energy tail in er rcg (cj) seen in figure 7(a) for all n is due to 
incoherent excitations of type D (figure 1), which have a maximum energy of 2(W + uj ) 
(within the WC approximation). 

Strong coupling Figure 2 illustrates the contributions to a rcs (uj) in the SC regime, where 
the polaron band lies entirely in the interval (—uj ,uj ). Equations (A1)-(A3) represent 
transitions between the states in the polaron band accompanied by multi-phonon processes. 
The processes (B) and (C), corresponding to transitions between the incoherent polaron 
spectrum and the polaron band states, again involving multiple phonons, are expected to 
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Figure 9. Analytical results for a rc s(ui) ( ) and 5 rog (o;) ( ) from the SC 

approximation. Here u/o/t = 4, .Ep/t = 8 and n = 0.4. 

become important for intermediate coupling. 

Examining the analytical expressions, the most important processes in the true SC limit 

2 

(A, g > 1) appear to be those of type (A3), with the dominating (with respect to e~ 9 ) 
contribution to a ICS (ui) given by equation (B.14). In this case, o- Icg (uj) consists of sharp 
peaks at iv = suj with heights modulated by the Poisson distribution with parameter 2g 2 , as 
may be seen in figure 5(c). The maximum of the conductivity lies near the maximum of the 
Poisson distribution determined by s = 2g 2 = 2E p /uj Q . 

The peak heights in the anti-adiabatic SC case (g 2 — 2) shown in figure 9 can 
be understood by evaluating the relative weights in front of the integral on the rhs of 
equation (B.14), yielding 1, 1, 8/9 for s = 1, 2, 3, respectively. 

The picture of the optical conductivity implied by our findings agrees with results 
obtained previously in the framework of small-polaron theory at T = [13,42]. The 
separated peaks may be smeared out by means of an additional "smoothing mechanism" [42] 
or due to limited experimental resolution. The evolution of the small-polaron absorption 
spectrum towards a smooth function due to damping and finite temperature effects was 
considered in [43], and Emin [39] has presented the absorption spectrum in the form of a 
Gaussian centered at 2E p , assuming the broadening of the energy levels of the localised 
polaron states. 

Crossover from weak to strong coupling Figure 10 illustrates the crossover from weak to 
intermediate to strong coupling for a fixed density n = 0.4 in the adiabatic regime ui /t = 0.4 
within the IC/SC approach. The corresponding spectral functions can be found in figure 14 
of [17], and we use the same values of the parameter R which decreases with increasing E p 
reflecting a decrease of the polaron size. 

The position of the maximum in cr rog (a') changes significantly as a function of E p . For 
weak coupling, it is located just above the absorption threshold at ui = ui Q (figure 10(a)), 
whereas for E p /t = 2 it lies near 1.5u Q . Nevertheless, the frequency range for absorption 
in panels (a) and (b) is very similar. For E p /t = 4, similar to the low-density case shown 
in figure 5(c), er rcg (o>) peaks near the small-polaron value 2E p . The analytical approach 
overestimates the band narrowing, yielding SC behaviour even for a value E p /t — 4 lying at 
the border to the SC region. Since for this value of E p the IC and SC approximations yield 
virtually identical results, we have used the SC approximation to obtain a consistent value of 
grcg -p ne s i m ii ar ity of the SC results for n = 0.1 in figure 5(c) and n = 0.4 in figure 10(c) is 
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Figure 10. Analytical results for a rc s(u)) ( ) and S rcg (w) ( ) from the IC [(a),(b)] 

respectively SC [(c)] approximation. Here u>o/t = 0.4, n = 0.4, and (a) Ep/t = 0.1 
(R = 1.5), (b) E P /t = 2(R = 1.3) and (c) B P /t = 4. 

due to the rather weak density dependence of our analytical results pointed out before in the 
case of the spectral function [17]. 

The evolution of S rcs with increasing E p is similar to that observed in the one-electron 
case (figure 3(a)), with a maximum in the IC regime. Whereas for weak coupling transport is 
mainly coherent (as reflected by a large Drude weight £D/<j q w 0.61 ^> tS res /a ), S rcs 
increases noticeably for E p /t = 2 owing to enhanced incoherent scattering. For strong 
coupling, since V — > 0, the reduced kinetic energy suppresses S rcg via the f-sum rule. 
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Figure 11. Analytical results for a Ic s(ui) ( ) and S Ies (uj) ( ) from the IC 

approximation. Here uj /t = 0.4, E P /t = 2.0, (a) n = 0.1 (R = 2.1) and (b) n = 0.3 
(R= 1.5). 



Density dependence at intermediate coupling We finally examine the density dependence of 
the optical response in the IC regime. To make a connection with previous work [16], we take 
uj Q /t — 0.4 and E p /t = 2. Figure 11 depicts results for a rcg (uj) for n = 0.1 and n = 0.3, 
whereas the case n = 0.4 is reported in figure 10(b). 

The low-density case n = 0.1 has already been shown in figure 5(b), and results are 
presented without rescaling in figure 1 1(a) for reference. Compared to n = 0.3 (figure 1 1(b)), 
we observe a transfer of spectral weight from high to low frequencies, causing a strong 
reduction of the broad high-energy hump. This trend continues upon increasing n even further, 
as can be seen from figure 10(b). Besides, the continuous increase of S ree with increasing n 
is similar to the WC case depicted in figure 4(a). 

These features of the optical conductivity suggest that the increase of the kinetic energy 
with n, as shown in table 1, has a more profound origin than the mere increase of the charge 
carrier density. In this connection, it is necessary to discuss the dependence of the "optimal" 
polaron radius R on density n. The values reported in the caption of figures 10 and 1 1 illustrate 
that within our variational approach, R decreases with increasing carrier density, i.e., R = 2.1 
for n = 0.1, R = 1.5 for n = 0.3, and R = 1.3 for n = 0.4. This seems to be at odds with the 
crossover from polaronic behaviour at small n to metal-like behaviour at larger n observed 
numerically [7, 8], since a decrease of R with increasing E p at fixed n generally enhances the 
polaronic character of the spectra [17]. 

Whereas the dependence of e kin on the coupling strength at fixed n may be interpreted 
as the crossover between the small-polaron and large-polaron regime, the concentration 
dependence of e kin seems to be more involved [7,8]. Indeed, according to the variational 
treatment of the IC case, the concentration dependence of the kinetic energy and optical 
conductivity is to be explained by the n-dependence of the parameters 7,7,77 defined by 
equations (5), (6). These parameters are determined, using the variational variables R, by the 
energy balance of the Hartree energy of polarons and the energy of the lattice deformation 
background [17], which at A = 1 may be quite significant. Although the optimal R increases 
with decreasing n, this energy balance leads to an increase of kinetic energy with increasing 
n, as the resulting increase of the lattice deformation is sufficient to enhance the mobility of 
charge carriers. 

This picture may be supplemented and improved by a more detailed insight provided 
by the electron-lattice and electron-electron correlations calculated numerically in [8]. 
According to these results, the increase of the carrier concentration is accompanied by the 
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short-range development of the charge-density wave, which is connected with the spreading 
of the inhomogeneous lattice deformation about the centre of the polaron. Consequently, the 
analytical calculations based on the Hartree energy — taking into account only the averaged 
energy of the lattice deformation background — cannot describe these features revealed by 
numerical calculations. Nevertheless, the optical spectra turn out to reproduce the main 
features found numerically [8, 16], and the dependence of the total energy on n is in good 
agreement with exact data [8]. 

5. Discussion and conclusions 

The results of section 4.2 suggest that the analytical approach captures the main features of 
polaronic systems with finite carrier density. Therefore, it is interesting to relate our findings 
to the polaronic excitations observed in manganites, for which no general theory is available. 
To be more specific, we consider results for thin films of La 2/ , 3 Sr 1//3 Mn0 3 (LSMO) and 
La 2 / 3 Ca 1 / 3 Mn0 3 (LCMO) containing a characteristic mid-infrared peak usually assigned to 
polarons [5]. The shape and position of this peak differs substantially for the two materials, 
and has been interpreted in terms of large and small polarons, respectively [5, 44]. Apart from 
the shift of the maximum in a ves (u) to higher energies when going from LSMO to LCMO, 
the peak is highly asymmetric in LSMO — with a sharp onset just below the maximum and a 
high-energy tail — and much more symmetric in LCMO. 

In small-polaron theory, the maximum in a ree (uj) occurs at 2E p [40], so that experiments 
can in principle be used to determine E p . However, for small-polaron theory to be applicable, 
the resulting value has to be consistent with the high-temperature activation energy of the dc 
conductivity, i.e., 4E a = E p [40]. Experimentally, it turns out that this is usually not the 
case, indicating that small-polaron theory is invalid. This problem has been pointed out in [3- 
5], and a many-electron approach valid at intermediate coupling was proposed [3,27]. The 
basic assumptions are that perovskite LSMO falls into the WC regime, whereas LCMO is 
characterised by intermediate to strong coupling, in accordance with the Curie temperatures 
T c w 370K respectively T c « 270K [27], and that the polaronic character of quasiparticles 
is enhanced in thin films [3]. Hartinger et al [5] compare various theories to their data, 
concluding that none is able to satisfactorily explain the lineshape of the polaron peak. 

Based on the results of our simple analytical theory, we can make the following remarks. 
The change of shape from asymmetric to symmetric and the position of the maximum is 
observed when comparing the WC and IC/SC regimes. In particular, the WC approach 
features a sharp onset and a long high-energy tail, whereas the IC/SC results are reminiscent 
of the Poisson distribution of phonons. These main features can already be observed in our 
exact one-electron results in figure 5. Apart from setting the order of magnitude of a Tes (u>), 
the many-polaron effects discussed above are of minor importance for weak coupling (as 
appropriate for LSMO), but renormalise polaronic features toward "metallic" (or more WC- 
like) behaviour in the IC regime [7]. Although polaron effects are enhanced in thin films as 
used in [5], the discrepancy between E p as deduced from a ICS (ui) and E a for LCMO in [5] 
suggests that the sample lies near (T c « 240K) but not in the SC regime. The true SC case 
seems to be realised in the layered manganite La 1 2 Sr 1 8 Mn 2 7 , for which angular-resolved 
photoemission spectra [45] suggest the existence of small polarons even in the ferromagnetic 
state [3,46]. 

The present approach correctly predicts coherent states at the Fermi level giving rise to 
a Drude response, and reproduces standard small-polaron results in the SC limit. This is in 
contrast to the analytical treatment of the Holstein double-exchange model in [27], which 
includes the coupling of local and itinerant spins in the manganites. The WC theory for a gas 
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of polarons proposed in [19] may in principle explain the optical conductivity of LSMO, but 
it is not clear why long-range electron-phonon interaction is not screened in a dense metallic 
system. Exact results for the Frohlich model with one electron have been given in [47]. A 
model of optical absorption based on transitions between ground and excited mixed polaron 
states, formed by the hybridisation of the narrow small-polaron band, and the wide large- 
polaron band was presented in [42,48]. However, this model is dealing with a number of 
constants, not explicitly connected with the parameters of the Holstein Hamiltonian, which 
are to be determined by fitting the experimental absorption curves. Besides, it is not clear how 
the crossover from the SC to the WC case and the optical absorption in the WC regime can 
be described in this framework. 

To conclude, we have developed an analytical theory for the optical conductivity of finite- 
density polaronic systems, based on previous work on the single-particle spectral function. 
The spectra are surprisingly rich and in satisfactory agreement with exact results. Together 
with previous numerical results [16], the present work provides a better understanding of the 
optical conductivity of dense systems with large or small polarons. Finally, it can explain 
some of the generic polaronic features observed in the manganites, but our electron-phonon 
model is too simple to make quantitative predictions. 
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Appendix A. 

The fermion spectral function (of electrons or polarons) can be written as 

1 ImS(fc,w) 
^[^-(^+^)-ReS(fc,w)] 2 + [ImE(^)] 5 

with the fermion self-energy T,(k,uj), and the free band dispersion of electrons (polarons) 
£ fe = — Wcosk(£ k — — Wcosk). 

For the Hamiltonian (1), the coherent part of the spectrum, A c (k, u) — characterised by 
Im E(fc, u) = — is non-zero for |w| < u . In this frequency interval, 

A c (k,w) = z k S[w-(E k + V )], (A.2) 

where E k is the solution of the equation 

£ fc = £ fc + Re E(fc, (A.3) 

and 



zl x = 



l-[^Re£(fc,c)U £fc+I) 



(A.4) 



Accordingly, the incoherent part of the spectrum, denoted as A lc (k,u), is non-zero for 
\u\ > uj and can be calculated according to equation (A.l). 

The optical conductivity is determined by the electronic spectral function A e (k,u). 
However, in the SC case, the correct fermionic quasiparticles are small polarons, and A c (k, u) 
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has to be expressed in terms of the polaronic spectral function A p (k, u) via [17] 

A e (k,w) = e~^A p {k,oj) + e"^ ]T ^ 

«>i 

x ^ j^4 p (fc', u) — slu q )9(lu — slu q ) + A p (k', lu + suj )8(—lu — slu q ) 

k' 

For the numerical evaluation of equation (21), we define 

A< (k, u < 0) = e^ 2 z k 5[u; - (E k + r,)]6(w + u) + e"^ ]T 



(A.5) 



ill 
s\ N 



z k ,S[uj + suj - (E k , + r?)]0(-w - suj )9[lu + (s + l)u ] (A.6) 



x 
k' 

and 

A>(k, u > 0) = e'^z k S[cj (E k + v )]e{u> - w) + e"^ £^ 

S>1 



^ z k ,S[u - suj - (E k , + Tf)]0(u - siv Q )9[(s + l)w - uj] , (A.7) 



x 

k' 

as well as 

A[ c (k,uj ^ 0) = c-° 2 A^(k,Lu)8(±Lj - cj )+c-° 2 " 



S>1 



ill 

s! N 
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k' 

Equations (A.5)-(A.8) for the SC case are formally the same in the IC case, with g replaced 
by 75 throughout. 

The contributions to the integral in equation (23) are then given by 

,o 

A= duj'e(oj + oj')A<(k,w')A>(k,w' + u), 

J — OO 

,o 

B= du>'6(w + u/)Af(k,u/)A£(k,u/ + 

J —OO 

,o 

C= &Je{u + J)A^(k,J)A>(k,uj' + 

J —OO 

r° 

D= du/0(w + u/)A£(k,u/)A*(k,w' +u). (A.9) 

J — OO 
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Re cr{ui) = - - ; , , (B.l) 



The Drude formula for the low-frequency conductivity, 

V T 

7T 1 + (W) 5 

gives at w = for the dc conductivity 

ct(0) = Pt/tt . (B.2) 
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Here r denotes the relaxation time. In the limit 1/r — > 0, we obtain 

Reo-(w) = D5(w). (B.3) 

The standard way to deduce the conductivity at ui = using the Kubo formula is to take 
the limit u> — ► of Rea(ui) calculated for uj > [40]. However, equation (21) yields zero 
for < u < w . To derive the Drude-like singularity (B.3) for the formulation of the f-sum 
rule, we shall start with a low but non-zero temperature T. Moreover, the spectral function in 
equation (21) is assumed to have a finite width — caused by a very weak additional scattering 
mechanism — leading to 

A(k, J) = ^___A___ , (B.4) 

^ I w - ( E k + v)r + a 2 

where A = (2t) _1 . In the WC approach, k = 1, t] = —fx, whereas in the IC/SC approach 
k = cxp(— g 2 j 2 ), and rj is given by equations (3) respectively (5). 

After substitution of equation (20) into equation (8), the limit ui — > gives 

1 f°° 

lim Rco-(o)) = cr(0) = -4cr — V(sinA:) 2 / Au'[A{k, w')] 2 d w ,/(w') . (B.5) 

N J-oo 

The subsequent limit T — > Ois taken using [d w ,/(w')] T=0 = —6(uj'), which after substitution 
of equation (B.4) into equation (B.5) yields 



lim cr(O) = ^^^y(sin/c) 2 z 2 — % — . (B.6) 

t^o v ' ir 2 N ^ y ' k \{E k + rj) 2 + A 2 ] 2 



Converting the sum over k into an integral over E, we obtain 

where the integration is over the energy range of the coherent part of the spectrum. However, 
considering A — > 0, the integrated function is non-zero only in a narrow interval around 
E k + r] = 0, i.e., for k w k F . Assuming the slowly varying part of the integrand to be 
constant, the integral reduces to 

/> Rr|rAf ^°"' (B - 8) 

Comparison with equation (B.2) then yields 

r 2 1 

t = Vt/tt . (B.9) 



CT (°) = ^( sinfc F) 2 .. , 

k—k-p 

According to equation (B.3), for 1/r — > 0, the Drude singularity takes the form 



d k E k 



[Rea{u)] D - V5(lu) = ^(sinfc F ) 2 

7T 



k 



d k E k 



S(uj). (B.10) 

k— fcp 



In particular, in the WC approximation with 7 = in equation (25) and k = 1, 

D = 2te 2 az /£ sin/c F . (B.ll) 

To check the above considerations, we briefly examine the free-electron case. In the latter, 
5 reg — > 0, and the kinetic energy per site is given as 

£ ki„ = Jj (- 2 *cosfc) = -^-2tsmk F . (B.12) 
|fe|<fe P 
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Since z k — 1, equation (34) is fulfilled. 

In the SC limit g 2 > 1, the leading term of V is 

VftsWe 2 aaiak F , (B.13) 

having the form of the Drude term of free carriers with half-bandwidth W. 

Comparing this estimate of V to the integrated regular part of the conductivity S* rcg , we 
find that the Drude part is negligible in the SC limit, so that the kinetic energy is determined by 
S reK according to the f-sum rule (34). To demonstrate this, we shall not use equation (24), but 
instead start with [Re<7(w)] A following from equations (A. 6), (A. 7), (A.9) and carry out the 
integration over ui before integrating over k. The dominating part of [Re a(ui)] A for estimating 
S TeK in the SC limit (corresponding to the contributions (A3) in figure 2) reads 



ReaM = Ijc- 2 * 2 £ ££J1 £ dk'z k ,6(-E k , „) 

x f dk"z k „d(E k „ +ri)S(w- slo q - E k „ + E y ) , (B.14) 
Jo 

where we assumed \E k \ < uj for all E k in the SC limit. The integration over u gives 



n) 

dk"z k „0(E k „ + rj)- 1 • (B.15) 

SUJ o - h k" + h k> 

To estimate S TCg , we set z k ,, z k „ = 1 and assume E k ,,E k „ <SC sw„. Accordingly, the rational 
function in equation (B.15) is replaced by l/sw and the integrals over k', k" simply give 
fc F and 7r — fc F , respectively. The Fermi wavevector corresponding to the band population at 
T = is k F = nir. The resulting estimate following from equation (B.15) is thus given as 

= ^n(l - «)Vc- 2 «i = ^n(l -n)/ 1 -) .(B.16) 

The brackets symbolise averaging with respect to the Poisson distribution with the 
parameter 2g 2 , and the concentration dependence of equation (B.16) is in agreement with 
standard small-polaron theory based on the atomic-limit approximation [40]. 

Comparison of equation (B.16) and equation (B.13) reveals that the Drude term in the 
f-sum rule (34) is negligible in the SC limit, so that e kin is proportional to S rcg , 



'kin - »(1 " ») ( " ) ■ (B.17) 

W \ S I 2g 2 

To make another check of the latter result, we compare equation (B.17) to the kinetic 
energy of one polaron at the bottom of the band, obtained from equation (A. 3) as t[d k E k ] k=0 
in [49]. Keeping only the leading term with respect to e~ 9 , setting z k = 1 and neglecting E k 
in comparison to w as above, we obtain 

W \ S I 2g* 

Assuming «<1, e kin = nE kin , in accordance with equation (B.17). 
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